ASYMPTOTIC FORMULA FOR THE SOLUTION OF THE 

STOKES PROBLEM WITH A SMALL PERTURBATION OF 

THE DOMAIN IN TWO AND THREE DIMENSIONS 

THI HONG CAM LUONG AND CHRISTIAN DAVEAU 

Abstract. In this paper we consider the resolvent Stokes problem in the 
case there is a small perturbation of the domain caused by a perturbed 
boundary. Firstly, we prove that the solution of Stokes problem is continuous 
due to this small perturbation. Secondly, we derive the first-order term 
in the displacement field perturbation that due to the deformation of the 
domain. It is worth emphasizing that even though only the first-order term 
is given, our method enables us to derive higher-order terms as well. The 
f—^ , derivation is rigorous and based on layer potential techniques. 
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CN ■ 1- Introduction 

^T ' Let ri C M (d = 2,3) be a bounded domain with boundary dQ of class C^. 

^r^ . We consider the Stokes resolvent system 

—Au + Vp = An in Q 

V-n = m n (1.1) 

K> ' \ u = g on do,. 

j^ I with u{x) is the velocity of the fluid, p{x) is the pressure, and g{x) is the 

boundary condition. 

Denote by v the outward unit normal to dO,, let dO,s be a (5— perturabation 
of dU as follows: 

In the two-dimensional case, let p € C^{dO), we consider the perturbed bound- 
ary 50^ be given by 

dils = {x : X = X + 6p{x)i'{x),x S dQ}; 
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and in the three-dimensional case, we consider the parallel surface 

di^s = {x : X = X + 6i'{x),x G 5J7}. 

Here we consider the parameter 5 > which tends to 0. We assume that 
g is an analytic function in a small neighborhood of dO,, and denote by gs 
the extended function of g onto d^s- Then denote by us the solution of the 
perturbation problem 

I —Aus + Vp5 = Xus in fls 

(1.2) 

In this paper we verify the continuity of the solution us with respect to 5 
and derive the asymptotic expansion of {us — u)\qo as 5 tends to 0, where Q° 
is any closed subset of il n fi^. The asymptotic expansion is derived due to 
the theory of layer potentials and Fredholm's alternative, and the properties 
of small perturbation of an interface. In connection with this, we refer to 
recent works in the context of interface problems [1], [2], [7], the continuity of 
the solution due to a small perturbation of an interface [7] and layer-potential 
theory for Stokes problem [3], [4], [5], [6], [8], [9]. In the work of H. Ammari, 
H. Kang, M. Lim, and H. Zribi [1] , they derived the asymptotic expansion of 
the boundary perturbations of steady-state voltage potentials, and now in our 
work, we derive the asymptotic expansion for the Stokes problem with Dirichlet 
boundary condition. However, by the same method, we can derive asymptotic 
formula for the Neumann problem as well. 

This paper is organized as follows. In the next section, we introduce some 
notations for small perturbations of an interface of class C^ , review some basic 
facts on the layer-potentials and give representation formulas. In Section 3, 
we verify the continuity of the solution with respect to 6, and in Section 4, we 
derive asymptotic expansion for the displacement field perturbation in term of 
6. 

2. Definitions and Preliminary results 
2.1. Small perturbation of an interface. 

2.1.1. In the case d = 2. Let a,6 G M, with a < b, and let X{t) : [a,b] -^ M^ 
be the arclength parametrization of 00, namely, X is a C^-function satisfying 
\X'{t) = l\ for aU tG [a,b] and 

dQ:={x = X{t),te [a,b]}. 

Then the outward unit normal to dO,, I'ix), is given by z^(x) = R^jLX'(t), where 

2 

R-n is the rotation by —tt/2; the tangential vector at x, T{x) = X'{t), and 

X'{t)J-X"{t). Set the curvature r(x) to be defined by 

X"{t) = T{x)u{x). 

We will sometimes use p{x) for p{X{t)) and p'{x) for the tangential derivative 
of p{x). 
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Then, x = X{t) = X(t)+5p{x)i'{x) = X{t)+5p{x)R^TiX'{t) is a parametriza- 

^ ^ 2 ^ 

tion of dVls- We denote by y{x) the outward unit normal to dVts at x. Then, it 

is proved in [1] that V(x) can be expanded uniformly as 

oo 

V{x) = ^ (5"z^"(3;), with x = x + 6p{x)u{x), x G dCls, x G dQ, (2.1) 

n=0 

where the vector- valued functions i'"(x) are uniformly bounded regardless of n. 
In particular, 

z/°(x) = u{x), iy^{x) = -p{x)T{x), x E dfl. 

Likewise, denote by das{x) the length element of dQs at x, which has an 
uniform expansion (see in [1]) 



das{x) = \X'{t)\dt = ^J{l-6T{t)p{t)Y+6'^p''^{t)dt = ^6''a"{x)da{x), 

n=0 

(2.2) 
with X G di^s, X G dQ, and u" are functions bounded regardless of n, with 

a^{x) = l,a^{x) = —t{x)p{x),x G dfl. 

Let x,y £ d^s^ that is x = x + 5p{x)i^{x), y = y + 5p{y)i'{y), then 

x-y = x-y + 6{p{x)u{x) - p{y)v{y)), 

and 

\x y\ — \x y\ ii-pzo . .^j ~\~o . .^j )• 

\x — j/p |x — yp 

Denote by 

p/ N {x-y,p{x)v{x) - p{y)v{y)) , . \p{x)v{x) - p{y)v{y)\^ 

E{x,y) := ■ r^ ,CT{x,y) ■= \ r^ • 

\x — y\^ \x — y\^ 

Since 90 is of class C^, we can see that 

\E{x,y)\ + |G(x,y)|2 < C||X||c2(an)l|p|lci(9n), for all x,y £ dn, (2.3) 
and hence 

oo 

\x-y\ = \x - 2/1 \/l + 2SE{x, y) + 5'^G{x, y) = \x - y\ ^(5"L„(x,y), 

n=0 

where the series converges absolutely and uniformly. In particular, we can see 
that 

Lo{x,y) = l,Li{x,y) = E{x,y). 
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2.1.2. In the case d = 3. Denote by dasix) the surface element of d^s at x. 
Thanks to the results in [2] we get the following uniform expansion for das{x) 

oo 

das{x) = ^ 5"-a"'{x)da{x), with x = x + 6i^{x), x G dQs, x G OO. (2.4) 

n=0 

where o"" are functions bounded regardless of n, with 

CT°(x) = l,a^{x) = -2H{x),a^{x) = K{x), x G dQ, 

with H and K denote the mean and Gaussian curvature of dQ respectively. 
In the other hand, due to the parallel property of dil and dQs we have 

i^{x) = I'ix), with X = X + 6iy{x), x G dCls, x G di^. (2.5) 

Let x,y £ di^s, then x — y = x — y + 5{i'{x) — viy))-, and 

|£-r = |x-,l^(i + 2j ^"-^;"^"\7^^^^ +^^ '"(f^-l^^'' ). 

\x — y\"' \x — y\^ 

Denote by 

E{x,y) := (^-^;K^)-Ky)) ^g(,^^) ,= H-) - <y)\\ 
\x — yp \x — 2/1^ 

Since dO. is of class C^, we can see that there exists a constant C depending 
only on dVl such that 

\E{x,y)\ + |G(x,y)|5 < C, for all x,y G 9^^, (2.6) 

and hence 

oo 

\x-y\ = \x -yyi + 26E{x, y) + 6'^G{x, y) = |x - y| ^(5"L„(x,y), 

n=0 

where the series converges absolutely and uniformly. In particular, we can see 
that 

LQ{x,y) = l,Li{x,y) = E{x,y). 

2.2. The potential theory for Stokes resolvent system. We start to re- 
view some basic facts in the theory of layer potentials. We consider here 
A G C \ {z G C : Rez < 0,Im2; = 0}. The fundamental tensors {T,F) of 
the Stokes resolvent system (jl.ip can be obtained by the Fourier transform 
method in the following forms (see [3] for d = 2 and [4], [8] for d = 3): 
In two dimensions, we have 

(Vijix) =-^{5,^.ei(VA|x|) + ^e2(VA|x|)}, ^^ ^^ 



for Vx G M"', X / with 



ei(K) =Ko{k) + k ^Ki{k) 



K 



-2 



e2(K) = -Kq{k) - 2k~^Ki{h) + 2K"^ 

where Kn (n G No) denotes the modified Bessel function of order n, and 5ij is 
the Kronecker' symbol. 
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In three dimensions, we have 

r r,,(x) =-^{5,^.ei(-VA|x|) + ^e2(-VA|x|)}, 

for Vx G M^,x / with 



(2.8) 



^i(^) = E|^^" = -^(^)(i-^~^ + ^"') 



6-2 



n=0 



oo 1 2 

^2(e) = E t;^^" = «^p(^)(-i + 3^"' - 3^"') + 3e-'. 

n=0 

By using the notations x = x — y = (xi, ■ ■ ■ ,Xd) and r = |3;|, we introduce 
the stress tensor 5 associated to the fundamental tensors (F, F) and having the 
fohowing components (see in [8]): 

S,^,ix,y):=-Fj{mk + ^^^ + ^^^, i,j,k = l,...,d. (2.9) 

Note that we used the Einstein convention for the summation notation omitting 
the summation sign for the indices appearing twice. We will continue using 
this convention throughout this paper. Taking into account the relation (|2.8p 
and (|2.9p we obtain the following explicit forms (see in [8]) as follows. 
In two dimensions, for \/x, y G R^, x ^ y we have: 

Sijk{x,y) = ~:^ \6,k^dl{VXr) + {6k,^ - S,,^)d2{VXr)\ (2.10) 

_ 1 r x^ ^^^^ ^ 2d2{V\r) - 2) 

with 

di (k) = 2^2 (k) + 1 - 4k-2 , (i2 (k) = 2^2 (k) + Ki^i (k) - 4k-2 . 
In three dimensions, for Vx, y G M^, x ^ y, we have: 



Sijk{x,y) 



-i [s^k^d^i-VXr) - (5,,^ + 5,,^)d2(-^/Ar)} (2.11) 
-^ {^^(3 - 3di(-^/Ar) + 2d2i-VXr))\ , 



with 

'^2(^) = E t4w^" = ^-^(^)(^ - 3 + 6^-' - 6^-') + 66-2. 

n=2 V" + ^J- 

The pressure tensor A associated to the stress tensor S has the following 
components (see in [3] for d = 2 and [4], [5] for d = 3): 



-^i6ikXlnr-A^-^ + 8^) for d = 2 
-^(^(Ar2-2) + ^) ford = 3, 
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for \/x, y G M , x ^ y. 

Next we consider the single- and double-layer potentials associcated with stress 

and pressure tensors. Let dD be a closed surface of class C^ of a bounded 

domain D and (p = {(pi,--- ,<j)d) a vectorial continuous function on dD. For 

X G R \dD, we define the single-layer potential Vofp which has components as 

follows: 



:= / Tij{x-y)(l)j{y)da{y), 1 = 1,--- ,d; (2.13) 

JdD 

and the double-layer potential Wd0 has components as follows 



vhHx 

dD 



Wh(p{x):= -Sijk{x,yH{y)(t>j{y)da{y), i = l,--- ,d; (2.14) 

JdD 

where ^{y) is the outward unit normal vector to dD at the point y. 
Additionally, we consider the functions Qd4> and liofp as follows 



TiD(t>{x) := Jqjj Aik{x - y)vk{y)(t>i{y)da{y) 



(2.15) 



for x G 



The functions {Vd4>-, Qd<P) and (Wd'/', II/j^) are smooth functions in each of 
the domains M.'^\Dq and Dq, respectively, where Dq is the inner domain with 
the boundary dD. All these functions satisfy the following equations of Stokes 
resolvent problem: 

-AVdHx) + VQd^x) = 0, V • VdHx) = 0, , . 

-AWdHx) + Vni5<^(x) = 0, V • UD(pix) = 0, ^ ^ 

for X G R^dD. 

For our prupose, we also need to introcduce the principal value of the double- 
layer potential in a point xq of dD defined by the following formula 



Vj;)(t>{xo) ■.= p.v. -Sijk{xo,y)uk{y)(t>j{y)da{y) (2.17) 

JdD 

For (p G {C{dD))'^, the following trace relation for Wd holds (see in [6]): 

Wd(P\±={T:^I + ICd)(P a.e.ondn, (2.18) 

where yVD<P \- and Wd<P |+ denote the limits from inside D and outside D. 

From the previous section, the solution for the interior Dirichlet problem can 
be presented by the pure double-layer potential (see in [4]), then taking D = Q 
and Z:» = ^^5 in the fomulae (i2J3]l . (fZTD . (l2J5]l and (f2l7D . we get the bound- 
ary integral representation of the solutions for the initial Stokes problem (jl.ip 
and the perturbation problem ()1.2p are {Wncp, Uncp) and (Wn^f/', nn^^), respec- 
tively. 
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Additionally, from ()2.10p and ()2.1ip . we conclude that the components Sijk{x, y) 
Vk{y) of S{x,y)u{y) can be written in the following manner (see in [3], [4]). 

In two dimensions, for \/x,y £M?,x j^ y , we have 

S,M^,y)uk{y) = -^ \j^di{VX\r\) + (^ - 6/-^)d2{VX\r\)^ , (2.19) 



1 J ViVjr • u 
with 



2^ , ,.,4 {2d,{VJ\r\) + 2d2{VX\r\)-2) 



di (k) = 2^2 (k) + 1 - 4K-^ da (k) = 2^2 (k) + kKi (k) - 4k-^ 
And in three dimensions, for Vx, y € M^, x 7^ y, we have 

S.,,{x,y)My) = -^ {]^4(-VA|r|) - (^ + s/-^)d,i-^\r\)^ 

(2.20) 



1 j rirjr ■ v 



with 

,2 



i3-3di{-VX\r\) + 2d2{-VX\r\)) 



oiJ (J/ 2 1 \ 



(n + 2 
n{n^ - 1) 



n=2 

,2 



^2 (" + 2)' 

For further purpose, we denote the kernel —Sijk{x,y)i^k{y) of double-layer 
potential }Cd4> by Vij{x,y). We note that T>ij{x,y) behaves like 0{r~'^~^^) as 
r — )■ , which implies this kernel is weakly singular, so the integral operator is 
bounded on {C{dD)f (see in [11, p. 243] or in [3], [4], [6]). 

3. Continuity of the solution with respect to the perturbation 

In this section we will show that the solution us of the perturbation problem 
(|1.2p is continuous with respect to (5 as 5 tends to 0. And for the guarantee of 
the solvability and uniqueness of the solutions for (jl.ip and ()1.2p . we assume 
the following compatibility conditions 



g • nda = 0, and / gs • nsdas = 0. 
an Jaus 

From the previous section, the solution to the Stokes system (jl.ip can be rep- 
resented by the pure double-layer potential on dQ 



u'{x) = w^^{x)= -s,jkix,y)uk{y)My)d(T{y), xen (3.1) 
Jan 

with the density vectorial function (f) satisfies 

Wncp \^= i^I + ICn)<P a.e.ondn. (3.2) 
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Likewise, the solution to the Stokes perturbation problem ()1.2p can be rep- 
resented by the pure double-layer potential on dQs 



ul{x) = Wli^cpix) = -Sijk{x,y)uk{y)<pi{y)da{y), X e^s, (3.3) 

JdQs 
with the density vectorial function (/) satisfies 

Wns^\-={ll + ICns)^ a.e.ondns. (3.4) 

From the boundary condition we have Wncj) l-= 9 on. dQ and WugCJ) \-= gs 
on di^s- Consider us — u in any closed domain 0° C il n fi^, from the above 
formulae 

us{x) - u{x) = WnM^) - ^nH^), x G 0°. (3.5) 

Let ips{x) = X + 5p{x)u{x) and ips{x) = x + 5i'{x) be the diffeomorphisms 
from dil. to Oils in the case d = 2 and in the case d = 3 respectively. The 
following estimates hold. 

Lemma 3.1. There exists a constant C depending only on Q such that for any 
function 4> G {C{dQs)Y, we have 

||(/Cn,<^) oips- K,n{4>°'4>5)\\{c{dn)Y < C^\\(i>\\{c{dn,)Y- 
Proof. Fix X G (9il, we then have 

{V^^4>)o^l,s{x)-K.U4>oi^s){x) 



= p.v. / [-Sijk{x, y)vk{y)3s{y) + Sijk{x, y)T^k{y)](t>j ° i^&da{y) 
Jan 
here js denotes the Jacobian of ips, and has the approximation js{y) = l + 0{6). 

We can express Sijk{x, y) = Sijkix — y), then due to the mean value theorem 
we have the following expression. 
In the case d = 2, we have 

-Sijkix, y) = -Sijk{x,y)-5{p{x)u{x)-p{y)v{y))VSijk{x-y+6{p{x)v{x)-p{y)v{y))), 

with Q <e <5. 

In the case d = 3, we have 

-Sijk{x,y) = -Sijk{x,y)-5{v{x)-v{y))VSijk{x-y+e{u{x)-u{y))), < 6* < (5. 

Then combining with (|2.ip and (|2.5p we get 



- Sijkix, y)vk{y) = -Sijk{x,y)vk{y) + 6Tij{x,y) + 0{5), 
with 

Tij{x,y) = -{p{x)u{x) - p{y)v{y))V Sijk{x -y + 6{p{x)u{x) - p{y)v{y)) 

—Sijk{x — y)vk{y) in the case d = 2; 
Tij{x,y) = -{v{x) - v{y))SISijk{x -y + 6{v{x) - v{y))) in the case d = 3. 
Denote by 



Thi 



o ^s) = p.v. / Tij{x,y){(j)j o ^s)ix)da{y), 
Jan 
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we have the equahty 



(/Ch,</.) o Mx) = ICh{<^ o i;s){x) + 6T^{(t> o i;s) + 0{6'). (3.6) 

By using the remarks 12.31 and 12.61 as well as the decay behavior of Sijk{x,y) as 
|x — y| —7-0 (as mentioned in Section 2.2), we obtain that Tij{x,y) = 0{\x — 
y\~'^~^^) as \x — y\ ^ 0, which implies that Tij{x,y) is also weakly singular. So 
the integral T^{(f> o ips) is bounded on {C{dQ)y, and we have that 

WTnQ) ° '4^s)\\{c{dn))'i < C'll'^o V'<5||(c(5n))d, 
this completes the proof. D 

Lemma 3.2. There exists a constant C depending only on Vt and g such that 

\\4> ° ^5 - H{c(an)Y ^ ^C{g,n). 

Proof. Since the integral operator ^I + ^n is invertible on {C{dQ)) , combining 



with the trace relation (|3.2p . (|3.4p and Lemma 1 we obtain: 

\\4>°'4's-<l)\\(c{dQ)r < C\\{-I + fCn)Q) o il^s - <l))\\(c{dQ)r 

+C\\{ICns^) o^s- ^n{4> ° '^5)\\{c{dQ)r 

< C\\gs oil's- 9\\[c(dQ.)Y + C6\\4)\\^c{dns)Y- 
Due to the mean value theorem and the analyticity of (^ in a neighborhood of 
do,, we can easily see that there exists a constant C depending only on 17 and 
g such that \\gs o -0 - g\\(c{dn)Y - ^^(d, ^)- 
We also have that ||</'||(c(ar2i))'' ^ ll55ll(c(9C))'' because of the invertibility of 

Finally we get that \\(p o ^^ - (/'||(c(an))d < ^C{g, Q). D 

Theorem 3.3. Let il.° is any closed subset of QCi ^s, there exists a constant 
C depending only on g and $7 such that 

\\us{x) - n(x)||(c(no))d < 6C{g, 9). 

Proof. As the double-layer potential is continuous on 17°, by using the extreme 
value theorem, there is a point xq G f7° such that 

us{xq) - u{xq) = Wn,4>{xQ) - Wn0(xo) = max \\WnMx) - Wn0(x)|| 

Recall from Section 2 that we denote —Sijk{x,y)i'k{y) by Vij{x,y), we then 
have 

U5{xq) - u{xq) = / ['D{xo,y)js{y)^o'iPs{y) -T>{xo,y)(j){y)]da{y) 
Jan 



'^{xQ,y)[J8{y)-l](t)oi)5{y)da{y)+ I 'D{xo,y)[(f)0'il;s{y) - (f){y)]da{y) 
an Jan 

I [P(xo, y) - V{xo,y)]<l){y)da{y) := h + h + h- 
Jan 



+ 
Ian 
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It follows from Lemma 2 that 

l2<6C2{g,n). 
By abusing of the approximation js{y) = 1 + 0{5) we obtain 

h<5Ciig,n). 
And tdue to the analyticity of T>{x,y) we have that 

The proof is completed. D 

4. Derivation of the asymptotic expansion 

Firstly, we investigate the asymptotic behavior of IC^g(j) as 5 — ?• 0. Denote 
hy (ps = (p ° ip5 the vectorial function with components (psj (j = ^,d), by the 
change of variable, we can rewrite the integral ICngcf) as the following form 



K^hs4>{x) =P-v. -Sijk{x,y)i/k{y)<t)5,j{x)das{y). (4.1) 

J an 

Then using Taylor expansion for the kernel of double-layer potential S{x,y) at 

6 = when r = x — y ^ and combining with (|2.2p and ()2.4p . we obtain. 

In the case d = 2 

- Sijk{x,y)uk{y) = -Sijk{x,y)uk{y) + 

' V ' 

(4.2) 

E^" E E -((^(-)-(-)j^(^)-(^))" v5g.,(x,.)).,^(.), 

n=l m+p=n \a\=m 

^ V ' 

In the case d = 2> 

- Sijk{x,y)uk{y) = -Sijk{x,y)vk{y) + 

' V ' 

(4.3) 
E^" E E ~^^"^"\,"^'^^'^ '^SUx,y)H{y). 

n=l m+p=n \a\=m 

^ V ' 

■■=Krjix,y) 
Then we introduce a sequence of integral operator {IC^-(j)s{x))n£N defined for 
any ^s G (C(50))'^ by 

IC^'Mx)=P-V- E / K:;;ix,y)a'^iy)cPsAy)My), ri<0. (4.4) 

Note that /Cq = /Cn • By using the same arguments as in the proof for Lemma 1 , 
we have that the operators JCqC^s are bounded on (C((9il))''. 
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Theorem 4.1. Let A^ G N. There exists a constant C depending on N and O 
such that for any (p G {€{80,$))'^! 

N 
n=2 

where (ps '■= 4>° '4'5- 

We now investigate the asymptotic behavior of ICns4' as 5 — >• 0. One can see 
from Theorem 2 that for each integer N, (ps satisfies 

1 ^ 

i^I + ICn + Yl ^"^n)h + 0(<5^+^) = g o ^-^ on d^. 

n=l 

Now using the Taylor development for g(x) = {gi(x))^^Yd at 5 = we obtain 
the following expansions. 
In the case d = 2 

9^i^) = E E ^V"5.(x)(Kx)Kx))" := J2 ^""Gnx) (4.5) 

n.=0 \a\=n ' n=Q 

note that, the first three terms are given by 

GO(x) =5^(x), G}{x) = p{x)HxyVg,{x), GJ{x) = \{p{x)i^{x))^V^g,{x){p{x)v{x))- 

where the superscript T denotes the transpose of a vector. 
In the case d = 3, we have: 

oo r-^ OO 

9^(5^) = E E ^V"5.(x)(K^))" := E ^"^"(^) (4.6) 

n=0|a|=n ' n=0 

with the first three terms are given by 

GO(x) = g,{x), G]{x) = u{x) • Vg,{x), GJ{x) = ^i,.{x)fv^g,{x){i.{x)). 
Therefore, we obtain the following integral equation to solve 

N oo 

71=1 n=0 

We look for the solution of (|4.7|) in the form of power series 

N 
n=l 

The equation ()4.7p then can be solved recursively in the following way: Define 

/ = (h + ICr'G^ = {h + ICy'g = <p (4.8) 

and for 1 < n < A^, 

1 / n— 1 \ 

" (4.9) 



{h + /C)-i f G" + ^ /C"-P0P j 
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We obtain the following Lemma. 

Lemma 4.2. Let N £N. There exists a constant C depending only on N and 
Jl such that 

N 

n=0 

where (jP are defined by the recursive relation fi4.9^ . 

Now we derive the asymptotic behavior oi us — u as 6 — t- 0. From the 
boundary integral representation of the solution in ()3.2p . after the change of 
variable we rewrite it as follows 

y^iiAx)= I V,,{x,y)4>sj{y)das{y) (4.10) 

Jan 

where cps '■= (p ° "^5^ ^^^ {4'5,j}j=Yd ^^^ components of cps- 

Then due to the Taylor expansion of Sijk{x,y) with respect to y, combining 
with (1211 and (1231). we obtain 



For the case d = 2 
- Sijk{x,y)ljk{y) = -Sijk{x,y)uk{y) 



(4.11) 



+ 2^"^ z^ Z^ -\ ^ySijk{x,y)Vk{y)- 

n=l m+p=n \a\=m 



-DU^W) 



T-3 



Likewise, in the case d = 3 



Sijk{x, y)^k{y) = -Sijk{x, y)vk{y) + X] '^" X] — T^y — ^y^ijkix, y)yk{y) ■ 

n=l \a\=n 



a! 



(4.12) 
Thanks to ([2:2]) . (plij) and (fil^ . the formula ([JTO]) can be rewritten as 
follows 

^hM{x)= f (v,,{x,y) + f26-V^j{x,y)] (4.13) 

Uiy) + E ^"-^"(y) + ^i^''^')] [ 1 + E '^"^"(y)) My) 

N 

A,(x,y)(/.,(y)da(y)+E'J" E / ^^'(^,y)'/''(y)^''(y)'i^(y)+o('5^+')- 
rrl _„ ,1~r .Jan 



^^ n=l m+fc+5=n-^^^ 
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Note that jQ^Vij{x,y)(t)j{y)da{y) = >Vq(/>(x). We then define, for n E N and 
for 2; G 0°, the terms n„ = {Un)^^Yd ^^ fohows 

<(^)= E / v^{^,y)<y\y)<t^){y)da{y). (4.i4) 

We obtain that 

AT 

us-u = Y^ 5"n„(x) + 0((5^+i), X G 9.°. (4.15) 

n=l 

The remainder 0{5 ^^) depends only on N and il. 

Let us compute the first order approximation of u^ — tt exphcitly. Note that 
(fp = (j) where (f) is defined by ()4.9p and 



4,^ = {h + lC^)-\G^ + lCy). (4.16) 

Therefore, in the two-dimensional case, ui = {u\)^^y2 takes the form 

u\ix)= [ Vlj{x,y)(f>,iy)daiy) - [ V,,{x,y)T{y)p{y)(t)j{y)da{y) (4.17) 
Jdo. Jan 

+ I Vij{x,y)4>]{y)da{y). 
Jan 

And in the three-dimensional case, ui = {u\)-^y3 takes the form 

u\(x)= [ V}j{x,y)cl),{y)da{y)-2 [ V,,{x,y)H{x)<l),{y)da{y) (4.18) 
joc Jan 

+ / Vij{x,y)4>]{y)da{y). 
Jan 
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